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^ Abstract 

^ We clarify the structure of thermal noise in AdS/CFT by studying the dynamics of an equi- 

^ librated heavy quark string. Using the Kruskal extension of the correspondence to generate the 

dynamics of the field theory on the Keldysh contour, we show that the motion of the string is de- 
I— I scribed by the classical equations of motion with a stochastic boundary condition on the stretched 

"5 horizon. The form of the stochastic boundary condition is consistent with the dissipation on this 

surface and is found by integrating out the fluctuations inside of the stretched horizon. Solving the 

^ equations of motion for the fluctuating string we determine the full frequency dependence of the 

'— ' random force on the boundary quark and show that it is consistent with the frequency dependent 

^T) dissipation. We show further that the stochastic motion reproduces the bulk to bulk two point 

^ functions of the Kruskal formalism. These turn out to be related to the usual retarded bulk to bulk 

00 

propagator by KMS relations. Finally we analyze the stochastic equations and give a bulk picture 
^ of the random boundary force as a flip-flopping trailing string solution. The basic formalism can 

be applied to the fluctuations of gravitons, dilatons, and other flelds. 
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I. INTRODUCTION 



The AdS/CFT correspondence has been used extensively to study the properties of 
strongly coupled gauge theories in a controllable setting [H [21 [3] . The correspondence has 
led to important insights into the nature of strongly coupled plasmas. In particular, Refs. 
[H El El 13) [H] established that in a large class of gauge theories with gravity duals the ratio 
between the shear viscosity and the entropy is 



This was important because it showed that there exist certain theories which realize the 
small transport time scales needed to explain the elliptic flow observed at RHIC [9l [TOl 
ITTl [T2| [T3l [Til [T5l [T6] . Since this work on shear viscosity many other transport properties 
of strongly coupled plasmas have been computed using the correspondence. Of particular 
relevance to this work is the computation of the heavy quark drag and diffusion coefficient 
in A/" = 4 Super Yang Mills (SYM) at large Nc and strong coupling |T7l [TSl 

Most of the time, thermal noise is neglected in AdS/CFT. This seems at odds with the 
fluctuation-dissipation theorem and leads to some seemingly incorrect results from the cor- 
respondence. For instance, it predicts the absence of long-time hydrodynamic tails |20] . 
zero drag on mesons [211 [22l [23] . and the lack of Brownian motion of a quark string in 
AdSs [17]. In many cases, the effect of thermal noise is suppressed either by large or 
large A, and therefore these inconsistencies with field theory intuition were rationalized as 
an artifact of these restrictive limits. Certain transport properties such as meson transport 
[21] and momentum diffusion [TH [23 [26] are intrinsically related to the fluctuations. These 
transport rates were computed using the correspondence by computing the drag and ap- 
pealing to the boundary fluctuation-dissipation theorem to determine the fluctuation rate. 
A notable exception to this rule is the calculation of the momentum broadening of a fast 
heavy quark [25] [26] which used the Kruskal formalism to compute this diffusion rate in an 
out of equilibrium setting where the fluctuation dissipation theorem does not apply [25] . 

The purpose of the present work is to overcome these difficulties by working through the 
simplest possible system which should exhibit drag and noise in AdS/CFT. This system 
is the Brownian motion of a heavy quark placed in the A/" = 4 SYM plasma. A quark in 
AdS/CFT is represented as an open string stretching from the horizon up to a probe brane. 
A schematic of the AdS geometry together with the heavy quark probe and the stretched 
horizon (see below) is shown in Fig. [l] The various problems with large Nc and large A are 
easily summarized by the simple fact that, according to the classical equations of motion, 
the string does not move in the absence of external force. Clearly this is not the gravity 
dual of a Brownian particle in equilibrium with plasma. 

One expects that the black brane should emit Hawking radiation inducing random motion 
on the string [27] . Hawking radiation is consistent with the fluctuation-dissipation theorem 
[2HI I2H1 [Sni [SB [321 [33j, and this consistency should lead to the correct picture of drag and 
diffusion of a heavy quark in the boundary theory. This idea has been sugested by several 
authors but was never clarified [511 [3S] • When analyzing dissipation in the context of black 
holes, the notion of a stretched horizon plays a central role [6l [361 [31] and this surface is 
shown in Fig. [T] at r = r/i = 1 + e. It will turn out that integrating out the fluctuations 
within the stretched horizon will yield a stochastic equation of motion with the required 
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FIG. 1: A schematic of a classical string in AdSs corresponding to a heavy quark. The horizon is 
at r = 1 in the coordinates of this work. The stretched horizon is at r/^ = 1 + e and the endpoint 
of the string is at with ^ 1. Gravity pulls downward in this figure. 



noise at the horizon. 

To keep this work self contained we have provided a fairly extensive review of the necessary 
ingredients. Section |Tl] reviews the contour formalism of thermal field theory and shows 

reviews 



how Brownian motion arises naturally from this algebraic structure. Section IV A 



AdS/CFT in the Kruskal plane and generalizes the Kruskal/Keldysh correspondence to 



o" = and the ra basis which is more natural for the classical AdS/CFT setup. Section VA 



shows that bulk to bulk contour correlation functions respect KMS relations. This result 



which, while known [29| [30] , was not entirely appreciated. Finally, Section VB presents 
a derivation of the stochastic force on the horizon and the subsequent section analyzes the 
results. A readable summary of the results and the bulk picture of Brownian motion is 



presented in Section VI 



II. LANGEVIN DYNAMICS 

A. Notation and Overdamped Motion 

In this section we will review briefly the Langevin process. Consider a heavy particle 
moving with velocity x, subjected to drag —rjx and random noise ^ 

Mkin|^ + = e , mm) = 2Tr^ s{t - f) . (2.1) 

Here and below we will only write the equation of motion for the x component of the motion, 
= ^a;. We also have anticipated that the strength of the noise {^it)^{t')) is related to the 
drag 1] through the fluctuation dissipation theorem. is the quasi-particle mass including 
in-medium modification of the mass. The diffusion coefficient can be written as 

D = -. (2.2) 

V 

Let us recall how this formula follows. Consider the ultimate long time limit. In this limit 
the — Mkincu^ term may be dropped since it is proportional to frequency squared and the 
motion is overdamped 



3 



t = — oo 



t = +00 



t = +00 — ia 



i/3 



FIG. 2: The Schwinger Keldysh contour. The fields labeled by "1" live on the upper time ordered 
axis, while the fields labeled by "2" live on the lower anti-time ordered axis. 



i.e. the dissipation exactly balances the force. Solving for the position of the quark, squaring, 
and averaging, we find the squared displacement 

{x\t)) = 2^t. (2.4) 
In the diffusion equation a Gaussian drop of dye spreads out as {x'^it)) = 2Dt leading to the 



identification in Eq. (2.2). This overdamped limit where the drag force exactly balances the 



velocity will be central in discussing the membrane paradigm and black holes. 



B. Langevin Dynamics From the Contour 

Next we review how the Langevin equations can be derived from the real time path 
integral [Ml [39] . The purpose here is not to track down every intricacy as there are reviews 
for this purpose [40j. In the recent literature we have found Refs. [111112] instructive. For a 
heavy particle coupled to an equilibrated bath of forces the real time partition function is 

[Dxi][Dx2] e'^'^^^^Q^^^ e-«/d<2M°i2 ! dtij^i(t^)x^(ti)-i J dt2J'2{t2)x2{t2)\ ^ (^2.5) 

/ bath 

where the path integration is along the Schwinger-Keldysh contour [HI |13l HH US]. The path 
integral over the "1" type coordinates (the upper line in Fig. [2]) represents the amphtude, 
while the path integral over the "2" type coordinates (the lower line in Fig. [2]) represents 
the conjugate amplitude. The path integral over the vertical pieces of the contour (which 
is not explicitly written in Eq. (2.5) but is implied by Fig. |2]) represents the average over 



the thermal density matrix e~^^ . The choice of a is arbitrary and reflects the fact that the 
thermal density matrix is stationary and therefore the average can be performed either at 
future or past infinity. An analogous dichotomy exists in the gravitational setup as discussed 



in Section |IV A[ We have not written the path integral for the bath which is also evaluated 
along the contour. 

When the particle is very heavy, the action is large and the motion is quasi-classical. The 
medium forces are small compared to the inertial terms and we can expand to second order, 
average over the bath, and re-exponentiate to find 
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Here the matrix of contour ordered correlation functions is 



Gii{t,t') -Gu{t,t') 

-G2l{t,t') G,2{t,t') 



(2.7) 



where for instance Gi2(t,t') = {Fi(t)F2{t')) is the average of the forces over the partition 

function of the bath. The relation to the operator formalism is the following: we define the 
time dependent operators 

:F,{t) = e'^':F{0)e-'^' , J^alt) = e+*^(*-'")^(0)e~'^(*-'") , (2.8) 

and then the correlation functions are 

tGn{t,t') = (^Tri{t):Pi{t')) , (2.9a) 

«Gi2(t,t') = ( J-2(t')-^i(t)) , (2.9b) 

iG2i{t.t') = (^p2{t) h{t')) , (2.9c) 

^G22{t,t') = Ifp2{t)p2{t')) . (2.9d) 



The KMS condition relates the different time orderings so there is really only one inde- 
pendent function which can be taken to be the retarded Green function 

iGR{t) = e{t)([Ht),Hm) ■ (2.10) 

\ / bath 

Using completeness and KMS relations it can be shown that 

iGn{to) = + iReGRitu) - (1 + 2n)ImGi?(cj) , (2.11a) 

iG'22(cu) = - iReGRiuj) - (1 + 2ri)ImG/j(tu) , (2.11b) 

iGi2{Lo) = - 2ne'^''lmGR{uj) , (2.11c) 

zG2i{uj) = - 2(1 + n)e~'^"ImG'K(o;) . (2. lid) 

Hidden in these relations is the inter-relation between the drag and noise. Almost all 
applications of real time thermal field theory have relied upon a rewritten version of the 
path integral known as the ra formalism (see, e.g., Refs. [HI HSj HHl HZ]). (Of particular 
note is the next to leading order computation of the heavy quark diffusion coefficient in 
weakly coupled A/" = 4 SYM which relied exclusively on the ra setup |16].) Here and below 
we take a = and the average over the initial density matrix is entirely in the past. Since 
in a quasi classical limit the amplitude is not very different from the conjugate amplitude 
we define the retarded (r) and advanced fields (a) for the particle and also for the forces 

Xr = ^ , Xa = Xx-X2, Tr = ^ , = ^1 " ^2 . (2.12) 

Xa should be considered a small parameter in the classical limit fI7j. Effecting this trans- 
formation we find 

Z = i \DX ][DX ]e~* ^ dtM^XaXr-f dtdt'Xa(t)iGu{t,t')Xr(t')-lxa{t)Gsyu,{t,t')Xa{t') ^2 13) 
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We have defined the propagators 



Gsym(^,0 = {rr{t)Tr{t')) = ^ {{T{t) , T{t')}) , (2.14a) 



iGRit,t') = (JS-W^a(^')) = ^(^) ([^(t),^(0)]) , (2.14b) 
iGA{t,t') = {Ta{t)Tr{t')) = -d{-t) ([^(t),^(0)]) , (2.14c) 
made use of the famihar contour relationships 

Gsym = \ [Gn + 6-22 + ^12 + ^21] , (2.15a) 

iGr = [Gn - 6-22 + 6-21 - 0,2] , (2.15b) 

iGa = ^ [Gn - G22 - G21 + G12] , (2.15c) 
= Gil + G22 - G12 - G21 , (2.15d) 
and also have used the reality relation between the advanced and retarded propagators 

GA{t) = Gni-t) , GAito) = G^uj) . (2.16) 
The KMS relations become simply the canonical form of the fluctuation dissipation theorem 

Gsymico) = -(1 + 2n)lmGR{uj) . (2.17) 
In Fourier space the path integral is simply 

J[DXr][DXa]exp(^-i J ^Xai-Uj)[-M^UJ^ + GR{uj)]Xriuj)^ g-i / x„(-a;)[G3y^(a;)K(c.) _ 

(2.18) 

After Fourier transforming the Gaussian by introducing a noise variable, 

e~^/i^^°(-'^)['^=ym('^)l^"H = / [D^]e+*/^''(~'^)^('^) e~^^=y™^ , (2-19) 



the partition function reads 

, _1 r du, a-^^Mu,) 

Z ^ [BXr][BXa][B^]e 2 J 2- GsymM 



X exp (^-z J "^Xai-co) [-M^LU^Xricu) + GR{uj)xr{uj) - i{uj)] ^ . (2.20) 

At this point one may integrate over Xa{—uj) yielding the path integral 

[Bxrm]e '^--M 5^ [-M^QU^Xriu) + GR{u)xr{u) - ^{u)] . (2.21) 

This equation means that the partition function is simply an average over the classical 
trajectories under the influence of a random colored force, 

[-M^^2 ^ ^ ^ (^(-^)e(^)) = aym(^) = -(1 + 2n) lmGR{u) . 

(2.22) 



In time this equation reads 

MQ^ + f dt' t') X{t') = at) , mm) = t') , (2.23) 



which is a generahzed Langevin equation with the drag and corresponding noise 
One method to implement such colored noise on the computer has been given in Ref. [50] . 
At small frequencies we can expand the retarded Green function to u"^ 

Gji{uj) = -itori - AMu^ , (2.24) 

and then the effective equation of motion is the original Langevin equation 

Mkin^ +v^ = m , mm) = 6{t - f) , (2.25) 

where we have defined the kinetic mass as Mkin(T) = Mq + AM. 



III. REVIEW OF TRAILING STRINGS 

Our purpose here is to collect some of the results from the AdS/CFT correspondence on 
the drag and diffusion of heavy quarks [T71 [THl [IS]- A canonical choice of coordinates for 
the metric of black hole AdSs (which we will denote with bars) is 

rf.^ = g[-/(&f)d^^ + ^x^]+;|^^ (3.1) 

with /(r) = 1 — Here h is the inverse horizon radius which is related to the Hawking 

temperature, h = l/vrTL^. We will use a different set of conventions defining r = bf, such 
that r is a measure of energy in units of vrT 

ds' = ircTfL' [-r'f{r)dt' + r'd^'] + . (3.2) 

We are considering the dynamics of very massive quark which is represented in AdSs by 
a long string stretching from the horizon upwards towards the AdS boundary terminating 
at Tm as illustrated in Fig. [T} This straight string is a solution to the classical Nambu-Goto 
equations of motion. Then we consider small fluctuations around this long straight string. 
The action for these fluctuations is (see Appendix |A]) 

(3.3) 



S=- I dtdr 

where 



m + ]^To{r){drxf - ^^{dtxf 



27r£2 

has the meaning of the local tension, and 



Ur) = fr' = i^/Xn'T'm fr' , (3.4) 
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is mass per unit r. The zero temperature mass of the quark Mq = mrm- The speed of waves 

on the string is Cs = \/To{r)f /m = irTfr"^ and therefore waves propagate from Vh = 1 + e 
to the boundary in a time of order ~ l/7rTlog(l/e). 

For a prescribed boundary value Xo{oo) we can solve for the classical waves on the string 
by imposing retarded boundary conditions at the horizon. The classical equations of motion 
after Fourier transforming in time are 

2 

dr{To{r)drX{uj, r)) + —x{uj, r) = . (3.6) 

We will define F^{r) as the retarded boundary to bulk propagator, i.e. the solution which 
satisfies 

lim F^{r) = 1 , 

and retarded boundary conditions. 

The waves on the string due to the sinusoidal motion of the end point of the string are 

Xo{uJ, r) = Xo{uj)F^{r) = Xo{uj) + [tan^^ (z) - tanh"^ {z)] + 0{uj^) , (3.7) 

with z = 1/r. The term multiplying the velocity —iuJXo{uj) is the "trailing string" solution 
of Refs. pm [T9] . (In writing this equation we have assumed that we are not exponentially 
close to the horizon, i.e. ci;log(l/e) -C ttT.) Thus to lowest order the string trails behind the 
sinusoidal motion according to the expected form. The next term in this series is proportional 
to the acceleration and is analyzed in Appendix |Xj 

The retarded force-force correlator is found by taking the boundary limit 

TO, ,2 



G"^ = lim T,(r)F_^9,F^ = -M^u'' + Gnico) , (3.^ 



where Mq is the zero temperature quark mass, and the term —MqUj"^ arises from the "diver- 
gent" part of the boundary limit. Then using To(r) given above and the retarded solution 



given in Eq. (3.7) and Appendix |A| we find 

GR{ijj) = -iuri - AMlo^ , (3.9) 

with 

1 



f] = -VXnT\ (3.10) 



AM = -^. (3.11) 
As discussed above, the field theory interpretation from this form for the retarded force-force 



correlator is that the quark obeys the stochastic motion in Eq. (2.23) with the specified 
transport coefficient r] (first computed in Ref. pTl [T8| [19] ) and in medium mass shift AM 
(first computed in Ref. [IZ]). The primary aim of this work is to show how this stochastic 
equation is derived in AdS/CFT and to give a bulk picture to the stochastic process. 

Before taking up this enterprise we make the following technical notes about the solutions 



to the bulk equations of motion Eq. (3.6) 
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1. The conjugate solution satisfies F* = and obeys advanced boundary conditions. 
When constructing solutions we note that ImF^(r) obeys the same differential equation 
but is normalizable, i.e. ImF(r) — > as r oo. 

2. The imaginary part of the retarded Green function, 

ImG^H = ^ [F_^(r)9,F^(r) - F^(r)a,F„^(r)] , (3.12) 

can be evaluated at any radius. This is because the term in square brackets is the 
Wronskian of the differential equation which depends on r in a particular way which 
precisely cancels the leading To(r) factor. This fact will be used repeatedly in Sec- 
tion |VS1 

3. When evaluating the evolution of perturbations in the bulk, it is often very useful to 
define the retarded bulk to bulk propagator. Many authors {e.g. Refs. [52l [53] ) have 
defined a Green function which is infalling at the horizon and normalizable at the 
boundary 

lmF^{r)F{f)9{r,f) + F^{r)lmF^{f)9{f,r) 

where we have defined the retarded Wronskian x To (r) 

To{r)W,etir) =To{r) [ImF (r)F(f) - F'{r)\mF{r)] , 

= - ImGniuj) . (3.14) 

Note again this combination is proportional to the retarded force-force correlator and 
is independent of r. 



IV. THE BOUNDARY PICTURE OF STOCHASTIC MOTION 



A. AdS/CFT on the Contour 

In this section we review how the Schwinger-Keldysh formalism is constructed in 
AdS/CFT by using the full Kruskal structure of the black hole. Here we will extend the 
results of Ref. [M] only slightly to show how the arbitrary a of the contour is present in the 
Kruskal formalism. We will then specialize to cr = and indicate how the results appear 
in the ra basis. The Kruskal/Keldysh correspondence was based on the pioneering work on 
Hawking radiation [2H1 EHl EHl EB E21 133] • 

The Kruskal plane of the eternal black hole is exhibited in Fig. |3] (see Appendix |B] for 
a summary of the Kruskal conventions adopted here.) Outside of the horizon there are 
two causally disconnected space time geometries both of which are asymptotically AdS. The 
right quadrant corresponds to the amplitude of the CFT while the left quadrant corresponds 
to the conjugate amplitude of the CFT. The dynamics of the two AdS/CFT's are coupled 
through boundary conditions at past and future infinity. 

AdS/CFT in practice amounts to a four step procedure. 
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FIG. 3: The full Kruskal plane. The right quadrant corresponds to the amplitude of the field 
theory (the "1" axis) while the left quadrant corresponds to the conjugate amplitude of the field 
theory (the "2" axis). 



1. First determine a solution to the classical equations of motion with certain boundary 
values. 

2. In order to uniquely specify this solution one must specify horizon boundary conditions. 

3. This solution is then substituted into the classical action which reduces to a boundary 
term. 

4. Finally, variation with respect to the boundary values determines the correlation func- 
tions of the CFT. 

In the full Kruskal plane these standard AdS/CFT steps are the following. 

1. First the equations of motion for the fluctuations on the string are solved in the 
full Kruskal plane with the boundary values 

lim x{uj,ri) = x°{u!), (4.1a) 
lim x{u!,r2) = X2{ui) ■ (4-lb) 

Here ri and r2 are the radial coordinates in the right and left quadrants respectively. The 
general solutions in the right and left quadrants are 

x{iu,n) = a{u;)FUn) + b{uj)F:{n), (4.2) 
x{u, ra) = c{iu)F^{r2) + d{u)F:{r2) . (4.3) 



Since we have only two boundary constraints given by Eq. (4.1), to uniquely specify the 
solution we must specify horizon boundary conditions. 
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2. The horizon boundary conditions specify the connection between the solution in the right 
quadrant and the solution in the left quadrant. We note that near the horizon the retarded 
and advanced solutions behave as 

e-'"*F^(ri) ~ e-''^^°s(^), (4.4) 
e-*'^*Fj;(ri) ~ e+^^^°s(-^). (4.5) 

The horizon conditions proposed in Ref. [5l] are that the solution should be analytic in the 
lower V plane and upper U plane. This physical choice is based on the intuition that if 
the right universe is to represent the time-ordered amplitude of the CFT, then the infalling 
mode (F^) should be positive energy (analytic in in lower V) while the outgoing mode {F*) 
should be negative energy (analytic in upper U). For instance in a free scalar theory the 
Green function at zero temperature is 

zGn{t,t') = ^e-^p(*-*')^(t,tO + ^e+*^p(*-*')^(t',t) , (4.6) 

and therefore a one type source Ji{t') will produce fields with positive energy wave at later 
time and induce negative energy wave at an earlier time. With this choice Herzog and Son 
analytically continue from the right quadrant (V > 0, U < 0) to the left quadrant (V^ < 0, 
U > 0) yielding from the solution in the right quadrant x(ti;,ri) the solution in the left 
quadrant 

x{u,r,) = aiu) e-^"^F^{r^) + h{uj) e^-^l'^ F:{r^) . (4.7) 

Now the values of a{uj) and h{ijj) can be determined from the boundary values x^^iu) and 
X2(co'). The resulting solution ultimately reproduces the contour correlation functions for 
the specific choice of cr = /3/2 [M] . 

Here we will generalize their work slightly by extending V \V\e~^^ and —U 

original work 6 = n. With this choice for the analytic continuation 
of U and V, the radius is fixed since {—U)V is changed ultimately by a factor e~^'^* (note 



Eq. (B4)). We will find that this choice for the analytic continuation reproduces the contour 
relations for arbitrary cr. Loosely speaking, the difference between these choices is a decision 
about whether to perform the thermal average at past infinity ^ = or future infinity 
U = 0. We do not have a sharper physical explanation for this choice at this point. With 
this analytic continuation the infalling and outgoing solutions behave as 

FMi) ^ e-'^'^F^(r2), (4.8) 
Kin) - e+-/^e--'^F:(r2), (4.9) 



where we have defined^ 



The solution in Eq. (4.2) when analytically continued to the left quadrant behaves as 



xiu, ra) = a{uj) e-"'^F^(r2) + b{uj) e^"^'^ e-^"" F^ir^) . (4.11) 



Of course this definition will later turn out to correspond to the usual field theory definition. 
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Now as before we can solve for the coefficients a{uj) and h{uj) in terms of the and X2{oj) 
yielding the result 

a{uj) = a;^(tu)(l + n(cj)) -a;^(cu)e'^'"n(tu), (4.12a) 

b{uj) = x°2{uj)e'^''n{uj) - x1{uj)n{uj) . (4.12b) 



The solution is now fully specified by its right and left quadrant solutions, Eq. (4.2) and 



Eq. (4.11). 



3. Now we can substitute the solution into the boundary action 
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bndry 



—Xi{-uj,ri)drXi{uj,ri) H — 



2n 



X2{-UJ,r2)drX2{uJ,r2) 



and determine the generating function 
1 r duj 



(4.13) 



iS] 



bndry 



-xli-uj) [+iReG^ - (1 + 2n)ImGy x°{uj] 
+ x°2i-uj) [-iReG^ - (1 + 2n)ImG^] x^(cj) 



X 



li-io) [-2ne+""ImG'^] x°i 



u) -x^i-uj 



Here we identified the retarded Green function with 

F.^{r)drF^ 



\F{uJ,r) 



) [-2(1 + n)e~'^'"ImG'^] x°(cj) . 

(4.14) 



(4.15) 



4. Now by taking variations with respect to x" and X2 one can reproduce the full contour 
correlation functions which obey the appropriate KMS relations in Eq. (2.11). 



B. Simple Results in the ra Setup 



Before closing this review of Ref. [54] let us show how these results work out in the ra 
basis. Specializing to a = and introducing the bulk r and a fields 



xi{uj,r) + X2{uj,r) 
Xr{uj, r) = , Xa{uj, r) = xi[uj, r) - X2[uj, r) 



(4.16) 



we find that the bulk fields are rather simply related to the boundary fields in the ra setup 



Xa{uJ,r) =F*{r)xl{uj) , 

Xr{iv, r) =F^{r) x°{uj) + i{l + 2n) \mF^{r)xl{uj) . 



The boundary action in the ra formalism after rewriting Eq. (4.13) is 



5, 



To{rm) 



bnry 



2ti 



[xa{-uj,r)drXr{uj,r) + Xr{-uj,r)drXa{uj,r)\ . 



(4.17a) 
(4.17b) 



(4.18) 
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Substituting the solutions given in Eq. (4.17) gives the desired result in the ra setup 



iS] 



bnry 



with G'sym(u;) 



2^ 



duj 
2^ 



Xli-Cj) [Gsy^(^)]<(^) 



(4.19) 



— (1 + 2n)lmG Ftiuj) . Of course this could be found simply by changing 



variables in Eq. (4.14), but it is hoped that the relatively simple causal relations in Eq. (4.17) 



will be useful in understanding non-thermal thermal fluctuations of AdS/CFT in an out of 
equilibrium setting [55j. 



C. Stochastic Motion on the Boundary 

Now we will show how these results lead directly to stochastic motion on the boundary. 
A formal path integral for the string is 

Z = J [Dx^Dxi] [Dx^Dx2]e^^i-^^% (4.20) 

where in the right quadrant of the Kruskal plane, Dx" is a temporal path for the string 
endpoint, Dxi is a bulk path integral for the body of the string, and the analogous "2" 
symbols are defined for the left quadrant. Slightly more explicitly 

[Dx?] = Y[ dx°^{t) , [Dxi] = Y[ dx{t, r) . (4.21) 



t,r 



Imagine integrating out the bulk coordinates to find an effective action for the string end- 
points. In a Gaussian approximation the integrals can be done. The result is a determinant 
times the exponential of the action evaluated with the classical solution which passes through 
the two endpoints and which obeys the appropriate horizon boundary conditions. This is of 
course simply the boundary action 

Scs = --^^ir^f ^xi{uj,ri)drXi{uj,ri) + '^"^^"'^ [ ^X2{-uj,r2)drX2{uj,r2) , (4.22) 
2 J^^ 2tt 2 J^^ 271 

which can be evaluated using the results of the previous section 

tS:^ = -ij ^<(-^) [-M^qu' + Gn{u)\ <(^) {G.yA^)] • 

(4.23) 

Now the partition function for the string endpoint is 

Z = j [Dx°D<] e^^°ff . (4.24) 
At this point we recognize the same partition function as discussed in the introduction, and 



following the discussion after Eq. (2.20) we conclude that the string endpoint obeys the 
expected equations of motion 

[-M°qU^ + Gr{u)\ x\u:) = r(^) , (r(-^)r(^)) = G,,^{u) = -(1 + 2n) ImGniu) . 

(4.25) 
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V. THE BULK PICTURE OF STOCHASTIC MOTION 



In the previous Section we integrated out all degrees of freedom of the string, except for 
the endpoints, to obtain the effective action for the endpoint coordinates. In this Section we 
will integrate out only the degrees of freedom inside a stretched horizon. In this way we will 
obtain an equation of motion for a string with a friction and a noise acting on the string at 
the stretched horizon. 



A. Bulk to Bulk Correlators in Contour AdS/CFT 



In this section we wish to compute the bulk to bulk correlators in the Kruskal formalism 
of AdS/CFT. For simplicity consider correlation functions of disturbances along an infinitely 
long string with fixed endpoints stretching from the horizon to the boundary. The action 



for these fluctuations has been given in Eq. (3.3). We will consider the generating function 



Z[Jl,J2 



PXi] [DX2] e^^^~^^^ e* dtidri Ji{ii,ri) a;i(ii,ri) ^-i f dt2dr2 J2{t2,r2) X2{t2,X2) j^^ 



where as before [Dxi] indicates a bulk path integral. We wish to compute all bulk two point 
functions 

(5.2) 
(5.3) 



ss' ' 



with 



1 



52 InZ 



%^ 6Js{t,r)6Js'{t,f) 
To this end we note that the classical equations of motion are 

,2 



Gil —G12 
—G21 G22 



dr{To{r)drx{uj,r)) H —x{uj,r) + J{uj) 



0, 



(5.4) 



where J{uj) is the external force per unit length on the string. Our task is to construct a 
Green function of this differential equation 



TfliiJ 

dr{To{r)drG{r,f)) + —G{r,f) 



S{r — r) 



(5.5) 



Given two solutions of the differential equation g^{r) and (7<(r) the Green function is easily 
constructed 

Ql^^ ^ 9>{r)g<{f)6{r, f) + g<{r)gy{f)e{f,r) 



where the Wronskian is 



To{f)W{f) 
W{f) = g'y{f)g<{f) - g'^{f)gy{f) . 



(5.7) 



In the present case we wish to construct the Green function in the full Kruskal plane. We 
should find two solutions which obey the boundary conditions. These boundary conditions 
are that g>{r) should be normalizable as r — 00 in the right quadrant of the Kruskal plane. 



14 



Similarly (7<(r) should be normalizable as r — oo in the left quadrant of the Kruskal plane. 
In formulas 

g>{T"i) = Im-^a;('"i) (right quadrant) , (5.8) 

and 

g<{r2) = lmF^{r2) (left quadrant) . (5.9) 

Using the boundary conditions of the real time AdS/CFT we extend across the horizon in 
the Kruskal plane. Writing ImF(ri) = {F{ri) — F*{ri))/2i and analytically continuing (?> 



according to Eq. (4.8) we have 



g>{r2) = ^ [e-'"^F^{r2) - e^'^ e-''^ F:{r2)\ (left quadrant) . (5.10) 
Similarly we realize that the analytic continuation of g<:{r2) to ri is 

g<{ri) = ^ [e+'^-F^(ri) - e--/^e'^-F:(ri)] (right quadrant) . (5.11) 

Now different correlators can be evaluated in a straightforward manner. Evaluating 



Gn{ri,ri) using Eq. (5.6), Eq. (5.8), and Eq. (5.11), and for simplicity taking ri > ri, we 
have 



W{fi) 

React(ri, fi) + «(1 + 2n)lmG,et{ri, n) . (5.12) 



Here n{uj) = 1/ {e^^'^ — 1) is the thermal occupancy factor and G^et (r, r) is given in Eq. (3.13). 
In deriving this result we note the intermediate results 

Win) =e^^'^{l - e-/^)lTyret(ri) , (5.13a) 

W{f2) =e--^(e+-/^ - l)iw^ret(r2) , (5.13b) 

and the much used identity n{—uj) = — (1 + n(c<j)). 

Continuing in this way, we summarize these results for the bulk to bulk correlators 

iGu{r,f) = + iReG'ret(r,f) - (l + 2n)Imact(r,f), (5.14a) 

iG22{r,f) = - zReGret(r,f) - (1 + 2n)lmG,^t{r,f) , (5.14b) 

iGi2 (r, f)=-2n e+'^'^Im G^et (r, f) , (5. 14c) 

iG2i{r,f) =-2(l + n)e"'^'"ImGrct(r,f). (5.14d) 

These are the familiar spectral and KMS relations between the retarded correlator and 
various contour correlation functions. It is reassuring that these thermal relations between 
the bulk to bulk correlators arise so easily in the Kruskal formalism although the result is 
not particularly new or surprising [221 ED] • 

To show that these correlation functions are indeed correctly normalized, one can shift 
the fields in the usual way 



Xs{uj,r) = 5xs{uj,r) - J dr [iG{iu,r,r)]^^, Js'{uj,r) , (5.15) 
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and determine generating functional 

Z[Ji, J2] = Z[0, 0] exp (^-^j ^drdf J,(-a;, r) r, f)],,, J,,(a;, f)^ . (5.16) 

B. Bulk Picture of Brownian Motion 

We will now develop a bulk picture of the Brownian motion. Rather than integrating out 
the entire bulk to determine an effective action for the boundary point, we will introduce a 
stretched horizon and integrate out only the fields which are inside the stretched horizon, 

rh = l + e. (5.17) 

The path integral for the fluctuations of the string is 

Z = J [DxIBxiBx^] [Da;° Dxa Dx^] [Da;^i D^] e'^'~'^'e*^^"*^' • (5-18) 

As before for the right quadrant for example, Dx° is the temporal path integral of the 
boundary endpoint of the string, Dxi indicates the bulk path integral of the string, Dx'^ 
denotes the temporal path integral of the string point on the stretched horizon, and 3x1 
the bulk path integral inside of the stretched horizon. 5*1 is the action integrated outside of 
the stretched horizon while 5*^ is the action integrated inside the stretched horizon. 

In a Gaussian approximation the integrals over the bulk coordinates inside the stretched 
horizon can be done. The result is a field independent determinant times the exponential 
of the action evaluated with the classical solution which passes through xj'(u;) and xl^^uj). 
This solution should also obey the contour boundary conditions given above. Substituting 
a classical solution into the action reduces to a boundary term at the horizon which gives 
an effective action for the horizon dynamics 

S^, = -^^ lx\{-u,r)d.x\{u:,r)'^ + ^^ [ f,{-uj,r) drxl{u;,r)^ . (5.19) 



Going through the same procedure outlined above we have the same structure for the 
effective horizon action 

^S^, = / ^ ^'(-^) ^ri^) - I J ^^'(-^) [-(1 + 2r^)ImG^(^)] a:„^^) , 

(5.20) 

but the retarded correlator is to be evaluated and normalized at the horizon. Using the 
near horizon behavior of the retarded solution F^{r) ~ 



Eq. (3.4) we find 



\F[u,r)\^ 



itUT]. (5.21) 
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The full partition function is now 



(5.22) 



Now in Eq. (|5.22|) we have the two halves of the Kruskal plane coupled by the effective 

/e will show that tf 
^write the action in 

Xi{r,uj) + X2{r,uj) 



action of the horizon. We will show that the effect of this coupling is to give rise to thermal 
noise in the bulk. We rewrite the action in the ra basis^ 



Xr[UJ,r) 

and the action in the bulk is 



Xa[uj,r) = xi[r,uj) - X2[r,uj) , 



iSi - 182 = -i I 7;— dr 



To{r)drXa{-uj, r) drXr{uj, r) 



muj'^Xriijj, r)xa{—0J, r) 

7 



We also follow what is now a standard procedure by introducing a horizon noise 
with the associated statistics 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



We now integrate by parts and obtain two boundary terms, one terminating at the radius 
of the string endpoint r^, and one terminating at the stretched horizon = 1 + e 



iSi - 182 + i8]^Q = - M W ^a(-^) [To{r,n)drXr{uj, r)] 

Jr 27r 



2^ 



x'^{-uj) [-To{rh)drXr{uj,r) - lur^x^iXu) - i'^uj)] 



- I I ^drxai-uj,r) 



-dr{To{r)drXr{uj,r)) - 



7 



(5.27) 



Finally integrating over a;°(— cj), Xa(— u;,r), and x^{—uj) in Eq. (5.22) we are left with a set 



of stochastic equations as in the simple example given in Section |11 B[ We will record these 
in the next section. 



C. Analysis and Discussion 

The preceding analysis leads to three equations 

1. The boundary endpoint of the string obeys the deterministic equation 

-T„(rja,a;,(cu,r) = 0. (5.28) 
This is simply the Neumann boundary condition for the free end of the string. 

^ The label r for retarded should not be confused with the radial coordinate. 
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2. The body of the string obeys the bulk equations of motion 



dr{To{r)drXr{uj,r)) H —Xr{uj,r) 



0. 



(5.29) 



3. Finally, the horizon endpoint obeys the stochastic equation of motion 

To{rh)drXr{u;,r) + t{u;) = -iujr]x^{uo) , (e^(-cu)e''(cu)) = 7]uo{l + 2n) , (5.30) 

with n{uj) = l/(exp(co'/T) — 1). Here we remind that To(r) = {y/Xn'^T^ /2)fr'^ is the local 
tension in the string. The meaning of this equation is that the motion of the horizon endpoint 
a;^(t) is overdamped. The resistance —rjx^ exactly balances the applied forces which in this 
case are the pulling due to the string outside the horizon To{rh)drXr and the random force 
^^{oj) which comes from integrating out modes behind the horizon. 

We now have determined a stochastic equation of motion for the endpoint of the string on 
the stretched horizon. These stochastic fluctuations are transmitted by the dynamics of the 
string to the boundary. We should show that the boundary endpoint of the string obeys the 



expected equation of motion in its most general form, Eq. (2.22). We should further show 



that the symmetrized two point functions obey the expected dynamics previously computed 
using the Kruskal extension of AdS / CFT 



(Ax,.(— c<j, r)Axr(co', r)) = — (1 + 2n)ImG].et(i-^, f) 

lmF,„{r)lm,„F{f) 
= (1 + 2n)- 



(5.31) 



First consider the average motion of the string. Averaging over the noise we find the aver- 
age string coordinates obey the usual equations of motion together with retarded boundary 
conditions 

ri [47rT(l - l/r)dr (x(cu, r)) + iu {x^)] = , (5.32) 

where we have used the definition for To(r). This equation together with the fact that (xr) 
is a solution says that the average obeys the retarded boundary conditions, i.e. behaves as 

(1 - 1/7.)-*W4tT_ 

Next consider the behavior near the boundary of AdS r — > r^. Here the equations of 
motion guarantee that the solution is a superposition of the non-normalizable mode and the 
normalizable mode which we can choose to be F^{r) and lm.F^^{r) respectively 



x{u,r) = Xo{uo)F^{r) + 



ImF^^(r) 



(5.33) 



We also have specified the boundary value of the non-normalizable mode Xo{oj) and have 
recognized that ^°{uj) must be a stochastic variable since the retarded solution reproduces 
the average motion. We have chosen to divide ImF^(r) by —\m.GR{uj) so that the ^°{uj) has 
the interpretation as the random force on the AdS boundary as we will show now. 



Plugging this functional form Eq. (5.33) into the Neumann boundary conditions 



To{r)drXr 



(5.34) 
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yields the expected equation of motion for the endpoint 

[-Mqu^ + Gr{uj)\ Xo{uj) = r(^) . (5.35) 
Our task now is to show that {E,°{—uj)^"{uj)) obeys the fluctuation-dissipation relation 

{i°{-uj)i\u)) = -(1 + 2n)lT^GR{u) . (5.36) 
Returning to the horizon we have the stochastic equation of motion 

To{rh)drXr{uj, r) + ^'^{uj) = -iujt]x^\uj) . (5.37) 



Substituting Eq. (5.33) into this equation we solve for i°{uj) in terms of i^iyj). Our task is 
simplified by first recognizing the origin of the first term on the right hand side 



iuri = Toirh) 



{rh)drF^{rh) 



(5.38) 



so that the equation for ^° resulting from Eq. (5.37) is 
r(^) 



-Toirn) [FUrh)drlmF^ - lmF{r^,)drF^{rh)] + F^{rh)e{uj) = . 



(5.39) 



The term in square brackets on the left hand side is the Wronskian of the two solutions and 
when multiplied by To(r) equals +lmGR{uj) and we find 



(5.40) 



Averaging according to the horizon statistics in Eq. (5.30), recognizing its origin Eq. (5.38), 
and again using the Wronskian relation Eq. (3.12), we determine the statistics of ^° 

{e{-uj)C{u)) = -(1 + 2n)lmGR{u) . (5.41) 

Thus random force in the boundary theory obeys the expected statistics of the fluctuation- 
dissipation theorem. 

Finally we would like to compute the bulk two point functions 



{Axr{—u!, r)Axr{uj, r)) 



(5.42) 



where Axr{uj,r) denotes the deviation from the behavior in the bulk due to the motion on 
the boundary Axr = Xr{uj, r) — x"{uj)F^{r). This is straightforward using the decomposition 
in Eq. (5.33), and the boundary statistics in Eq. (5.41); the result is 



{Axr{—uj, r)Axr{u!, r)) =(1 + 2n) 



lmF„,(r)lmF,„(r) 



-lmGR{uj) 
1 + 2n)lmG'i.ct(^^, r, f) . 



(5.43) 
(5.44) 



This result is naturally the same as computed previously using the contour correlation func- 
tion in AdS/CFT. It is also neatly consistent with the bulk fluctuation-dissipation theorem. 
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. _ ?X r = r/, = l + e 



FIG. 4: Balance of forces on the stretched horizon. The resistive force —rjx^ precisely balances 
the random force and the tension To leading to overdamped motion. 

VI. SUMMARY AND THE PHYSICAL PICTURE 
A. Summary 

In the previous section we have shown that a stochastic boundary condition emerges 
on the stretched horizon after integrating out the fluctuations inside this surface (see 



Eq. (5.30)) 



\3 



To{r,)drx{t,r,)+e{t)=vi\t) {em'it')) = G^^^{t - 1') . (6.1) 

Here i] is the drag of the horizon (or the late time drag of the quark), is the random force 
on the horizon endpoint, and C,^{uj) obeys the horizon fluctuation dissipation theorem 

G'.VM = +(1 + 2^V^- (6-2) 

This equation has a simple physical interpretation illustrated in Fig. |4j The motion of the 
horizon x^{t) is overdamped, i.e. the resistive force —rjx^ exactly balances the string force 
To{rh)drx{t,rh) and the random horizon force C,^. There is no transverse acceleration. 

This stochastic force on the horizon when transmitted to the boundary leads to an equa- 
tion of motion for the endpoint of the string 

Mq^ + I' Gn{t - t')x°(t') = i\t) , {Cm\t')) = aym(t, t') . (6.3) 



which is the expected generalized Langevin equation in Eq. (2.23). Here ^87111(1-^) obeys 
boundary the fluctuation dissipation relation 

Cy^co) = -(1 + 2n)lmGRiu) , (6.4) 

and Gr{u) is the usual retarded force-force propagator computed using AdS/CFT. The 
random force on the boundary ^°(t) is directly related to the random force on the horizon 



In this section the retarded r label (as in Xr) is understood. 
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Denoting F(t,r) the usual retarded boundary to bulk propagator, i.e. F^{rm) = 1 on 
the boundary and ~ (1 — i/7-)-«'^/4t7" t^q^^j- ^j^g horizon, the relation between the horizon 



and boundary stochastic forces is (see Eq. (5.40)) 



j dt'F{t-t',r,)e{t') 



(6.5) 



It takes a time of order ~ l/(7rT) log(l/e) for the noise from the horizon to reach the tip 
of the string, e should be considered small but not exponentially small so this time scale is 
really l/vrT. 

We also have a picture of the fluctuations in the bulk. The coordinate of the string in 



the bulk are given by two pieces (see Eq. (5.33)) 



x(t, r) 



(6.6) 



which reflect the retarded response to the boundary motion x°{t') and a deviation. The 



deviation Ax(t, r) is a random variable obeying the statistics (see Eq. (5.43)) 

(Ax(t, r) Ax(t', r')) = G^ymit - t', r, r') , (6.7) 

where G'sym(^ ~ t',r,r') is the symmetrized bulk to bulk correlator. This correlator was 
computed using the Kruskal extension of the AdS / CFT and is related to the imaginary part 
of the retarded bulk to bulk correlator according to a bulk fluctuation dissipation theorem 



Gsjm{uJ, r, r') = -(1 + 2n) ImG^ctiuJ, r, r') = (1 + 2n) 



ImF^(r)ImF^(r') 



-\v[vGr{u) 

The explicit form for the fluctuation amplitude in frequency space is 

. , . cot \^ IrnFjr) 

where is the boundary force. 



(6.8) 



-IraGniu 



(6.9) 



B. The Physical Picture 

Here we would like to consider the small frequency limit where an explicit analytic form 
for the retarded function F^^ is known and has a simple physical interpretation in terms of 
trailing strings. A quark in equilibrium moves quite slowly 



but it takes a long time tr for this heavy quark to randomize its velocity 

tr — • (6-11) 

rj 1 
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(This can be seen by examining the Langevin equations and neglecting the noise.) The 
distance the quark moves over this relaxation time xr is 



The dynamics that is observed depends on how the spatial Xobs and temporal resolution 
scales Tobs of the measurement compare to these scales xr and tr. 

First consider the time period over which quarks moves with nearly constant velocity v 

^<i-obs<rR. (6.13) 

For a quark moving slowly on the boundary with constant velocity v it will trail behind it 
a trailing string at least on average 

(x(t, r)) = Xoit) + t;AxTs(r) , (6.14) 

with 

AxTs(r) = — [iajr\z) - tanh-^(2;)] , (6.15) 

ZTTi 



and z = 1/r. Here we have used Eq. (5.33) for retarded response to the boundary motion 



and the explicit form of -FI;(r) at small frequency, Eq. (3.7). The term Axts is the "trailing 
string" solution of Ref. [T71 [19]. The distance between the head of the quark x°{t) and the 
average body of string is of order 

vAx^s ~ ~ ^i/4y ■ (6.16) 



This only gives the average behavior of the string. In general there is an additional 



random component which in the small frequency limit is white noise; using Eq. (6.9) we 
have 

Ax = x{t, r) - {x{t, r)) = — ^(t) Ax,,{r) , {Cm°{t')) = 2Tr]6{t - t') . (6.17) 

V 

Thus we see that around this average trailing string there is a stochastic ensemble of trailing 
strings which flip-flop around the head of the quark. This is illustrated above in Fig. [5] To 
estimate the amplitude of this stochastic process, let us imagine implementing this process 
on the computer where one would take time step At and then the width of the Gaussian 
process would be {^°^°) = 2Tri/ At. Taking the width At to be of order the memory time 
scale 1/T, we estimate that the string fluctuates around the average trailing string by an 
amount 

This is larger than the average deviation vAx-^s from the endpoint since it is not suppressed 
by I/a/tv^. Thus the average trailing string is nearly straight and the noise consists of a flip 



flopping trailing string solution. Notice the minus appearing in front of the Eq. (6.17). This 



is physically correct. When the random force on the boundary quark is positive, the string 
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^ < ^obs < 'TR 



h / 


1 








, - - ' average / 







•^obs 



A1/4T 



r = 1 



xr < a^obs 



r = 1 



FIG. 5: (a) The physical picture that emerges when observing the quark on relatively short time 
scales 1/T <C Tobs ^ tr- Here we show three subsequent time steps, ii,t2)i3; at each time step 
the string fluctuates to a new "trailing string" giving rise to a random force on the boundary. The 
average trailing string is perceived as a drag, (b) The physical picture that emerges on very long 
time and spatial scales. The horizon difl'uses and the string is brought along. 



is out in front of the quark. When the random force is negative, the string trails behind the 
quark. 

Now let us consider the ultimate long time limit where the time scales and spatial scales 
on which we are observing the quark are large compared to the relaxation time and 
relaxation length xr 

Tobs > "Tr , Xobs > VthTR Z — V'^m . (6.19) 

In this limit one can drop the influence of the bulk on the horizon dynamics, i.e. discard 



the TodrX in Eq. (6.1) since it is averaged over many different boundary velocities over the 
observation time Tobs ^ tr. The string in bulk is perfectly straight on the spatial scales we 
are considering. Thus the equation of motion obeyed by the horizon is 

— = —. 6.20 

T] dt ^ ^ 



This is the overdamped diffusive limit discussed in Section II A The result is that the horizon 



endpoint and the boundary endpoint diffuse in lockstep according to the expected rate 

2T 

{[x\t)r) = —t. (6.21) 



In summary, we have exhibited the full structure of the thermal noise on a fluctu- 
ating string in AdS/CFT. It is hoped that this will lay the groundwork to study the 
fluctuations of gravitons and other fields using the correspondence. The challenge now is 
to use the real time formalism in a truly out of equilibrium setting such as studied in Ref . [H] . 
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Note Added. In the final two weeks of preparation, a preprint appeared |56] which 
addressed some of the aspects studied here. In particular the stochastic boundary condition 



given in Eq. (6.1 ) was determined and the basic conclusion that the horizon diffuses according 



to Eq. (6.21) was also reached in this work. It is interesting that the derivation of the 



stochastic boundary condition in Ref. [56] differs substantially from the presentation given 
here. Presumably by marrying the two derivations a deeper understanding of Hawking 
radiation will be reached. 



APPENDIX A: SMALL FLUCTUATIONS AND THE TRAILING STRING 
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For small fluctuations we parameterize the string as 

(r, a) ^-^ (t = T,r = a,x = x(t, r)) 
determine the induced metric, 

\2 



'2 



The purpose of this appendix is to establish notation and to collect prior results. The 
Nambu-Goto action is 

S = -^J^ f drda \/-det hat ■ (Al) 

(A2) 

(A3a) 
(A3b) 
(A3c) 

(A4) 
(A5) 
(A6) 



hu = [-fr' + r'x']L\nT) 
1 

htr = [{x)x'r\nTY] L 



and write the action for small fluctuations as 



l + linT^fAxr-lj 



In the body of the text we write this as 



S=- dtdr 



771 

m + -To{r){drxf~—{dtxf 



where the local tension is 



and the mass per unit r 



m 



l-Kil 2 



(AT) 
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Then the equation of motion is 

-X + dr{fr%x) = . (A8) 



/ 

where, as is customary, we have defined tn = u/{ttT). 

This is a second order differential equation and there are two solutions. Near the horizon 
r — *• 1 the solutions are either infalling (-) or outgoing (+) 



1 



T A 



^--.) ■ m 

The solutions near the boundary r ^ consist of a normalizable and a non-normalizable 
mode. The retarded solution F^^{r) is the solution which approaches one near the boundary 
and is infalling at the horizon. More specifically near the boundary behaves as 

where the ellipses denote terms suppressed by additional powers of 1/r. In evaluating the 
motion of the quark we often evaluate the combination 

lim To{r)F.^drF^{r) = -M°uj^ + Gr{uj) , (All) 

where Gji{u) = (y^n'^T^ /2)B(uj) is the retarded force-force correlator. The first term MqUo"^ 
comes from the "divergent" rD^/2r^ term of the real part of the retarded Green function and 
we have identified 

as the zero-temperature mass of the quark where refers to the "canonical" coordinates 
). That this is the zero temperature mass can also be seen from the leading term 



in Eq. (3.1 



of Eq. (A4). 



In general Eq. (A8) can not be solved exactly. However we can set up a perturbation 
expansion at small frequencies, 

F^(r) = (1 - l/r4)-«'/4 [1 - zraF«(r) - w'F^^^r) + ...]. (A13) 

Substituting this ansatz into the equation of motion we end up with a hierarchy of differential 
equations. We solve these equations order by order in u by demanding that the solution 
behaves as 

F^{r) = (l -j X (regular function at the horizon) . (A14) 

We find 

Fi^\r) = - - In (^ + 1) - - In ^ 1) + - arctan (z) , (A15) 



25 



where here and below z = \jr. At quadratic order we have 

j.(2)(^) ^ J^- ^^ z{l-z)[z\n{z'' + l)-2z arctan (2:) + 2 In (z + 1) z - 4] ^ ^^^^^ 

This solution has a simple physical interpretation. The leading order solution is multiplied 
by the position of the tip string Xo{uj). Provided we are not exponentially close to the horizon 
we can expand the leading (1 — 2;^)"*™/^ factor. This yields the solution 

Xo{uj, r) = Xo{oj)F^{r) = Xo{uj) + v{uj)Axts + a{uj)Axa , (A17) 

where v{uj) = —iojXoiuj) is the velocity of the endpoint, and a{uj) = {—iuYxoiuj) is the 
acceleration of the endpoint. Further we have defined 

Axts(^) =77^ [tan"^(2) -tanh"^(z)] , (A18a) 



Axa(z) 



27rT 

^ ^ ^ ~' (Al8b) 



^log^(l-z^) + ^log(l-^^)F« + F(2) 



The leading term is given by Aa;Ts which is simply the "trailing string" solutions of 
Ref. [T71 [19]. Naturally to leading order in the frequency of the sinusoidal oscillation one 
simply recovers that the string trails behind the head of the quark with the expected form. 

The subleading term is described by the acceleration Ax a- A graph of this function is 
given in Fig. [6] and has a simple interpretation. Consider sinusoidal oscillations: when the 
head of the quark is moving forward and undergoing negative acceleration, the body of the 
string travels ahead of the trailing string solution due to inertia. Thus —Axa should be 
positive as r — > cxD reflecting the fact that the displacement is 180° out of phase with the 
acceleration. This inertial effect is indicated by the l/2r^ curve in Fig. |6} The dynamics close 
to the horizon stems from expanding out the leading (1 — z^)~^^^^ factor. The (1 — 2;^)"*"'/'^ 
behavior near the horizon has the interpretation that the string endpoint on the stretched 
horizon is overdamped. 

Using the solutions given above, we can expand these functions close to the boundary 

tt)^ \ (iw - ro^^ 



F.(r) = ^l + — + ...j + ^^-^(^l + ...j, (A19) 

and determine the retarded force-force correlator 

To{r)F^^drF^ = -M^oo^ + Gr{uj) , (A20) 
= -M^u^ -iuri-uj^AM . (A21) 
Here we have defined the transport coefficient 

V=^-^/XnT\ (A22) 

first computed in Ref. [TTl UHl [19] and the in-medium mass shift 

AM = -^, (A23) 

first computed in Ref. [1^. The fact that the mass shift is negative stems from the over- 
damped motion of the horizon endpoint. Also used in the text is the kinetic mass 

Mkin = M°+ AM . (A24) 
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-0.3 -0.2 -0.1 0.1 0.2 

-(jtT)2 AX3 

FIG. 6: The deviation of the string position from the traihng string solution (in units of 
\a{t)\ / (ttT)'^) during slow sinusoidal acceleration of the boundary endpoint. Here the accelera- 
tion is negative and the quark is moving forward. More specifically we are plotting —(7rT)'^Axa 
given in the text. 

APPENDIX B: NOTATION FOR THE KRUSKAL PLANE 

In this appendix we establish notation for the Kruskal variables used in the body of the 
text. We first define r^(r) 



rAr] 



1 



Then U and V are defined by the relations 

1 



The near horizon behaviors of = t + and = t — are 



(Bl) 

(B2a) 
(B2b) 



t + 



AtiT 
1 

— t — 

Also note that —UV is a simple function of r 



log(r - 1) , 
log(r — 1) . 



(B3a) 
(B3b) 



-U)V 



Ar*nT 



,2tan-l(r) 



r + 1 



(B4) 
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